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ABSTRACT 


The ratio cf geostrophic drag coefficients C/Cg under 
various stability conditions is obtained as a function of 
meenarcsonm Number, f Grape cocsurecrent Cr is defined in 
relatiorn™ te™ithe* mechanical mixim® Vergth, anmeewts ratic to 
uae actual arar=cocriicrent C is related to the Rivemardson 
number. “Byeusine Gbservational recent data on the normal- 
ized logarithmic wind shear of Monin-Obukhov one can obtain 
Cae a 

Thetwriter is deeply indebted to Dr. F. L. Martin 
(Professor of Meteorolory) for his suggesticn of the topic 
and his continued help throughout the investigation and 


Gteeine the preparation of this paper. 
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i. @entrocuction. 

Ume or the impempantt Crovlems 19 mlebecorolory is to 
esteblish relaticnships between the turbulent characteris- 
tics of the lower atmosphere and the large-scale synoptic 
paremeters. A significant turbulent characteristic is the 
@rcund Graz, which enables oné to specify the wind prorile, 
tive eddy Ciffusivity and the local energy cissipation in 
the surface layer, by knowing the topographical characteris- 
mrcs 1.,e€., tne romsnneds Penererer of & PATULCUIar =o lace, 

Lettau [7] had made an attempt to relate the drag 
coefficients with the static stability represented by a 
Gimensionless parameter called Richardson number. His 
investirations were based on 1955 O'Neill, Nebraska, data. 
[fireis the main purpose OF “tres a72oe. to Verily 2ane= exe 
mvs wnvestications using Igoe C’Neill data. It iS igereewe 
to note that this wrk deals with a relatively smooth 
surface, tierefore the "zero-plene displacement" which is 
sometimes used in the adiabatic-wind profile has been 
neglected, a procedure which greatly simplifies the com- 


putations. 





a. time Nature of the Data. 

The data used in this study were obtained by the 
Peeyect Prairie Crass [1] conducted meawee wel lly Nebraska 
during the summer 1956. Since this study rakes vse of some 
of the published data, i1.¢., wind speed, temperature and 
Sueeceifans Prt WS Mmeeessearnto destriWe Seme aspects of 
the instrumentation. The instrumertaticn employed was that 
of the mobile micrometeorolorical station of the Texas A&M 
preup. The station bas a slender aluminum mast supporting 
Six anememeters at he@ehts ofS, 4, 2,M1, 0.5 anae@veo meters, 
A similar mast supported seven temperature-mcasuring, radia- 
tion-shielded copper-constantan thermocouple junctions at 
ijemeints G,94,-2%8 2, Ord, 0725"angeO,125 meters. In addition 
the published report of Project Prairie Grass includes 
sectional sea level pressure maps at times nearly synoptic 
with the micrometeorolopical data. These maps revealed some 
small-scale features of the circulation near C'Neill which 
did not appear on the larger-scale facsimile maps. The micro- 
meteorolcrical data periods were selected to be nearly simul- 
taneous with the available maps. At times, linear interpola- 
tion between the available charts was used for computations 
OL =F ecestrephnic wind speeds. 

Regarding the errors of measurement, two’ types of teme- 
perature errorveccur: ise., callvratronmane racYaticon erircer, 
An estimate for calibration error was 0.0509C. The estimated 
radiation error was 0.1°C with the understanding tkat, in 
daytime with clear skies, and low windspeed, all measured 
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air temperatures overestimate the true temperature. At night 
with clear skies the measured temperatures wovld te lower 
than the ambient temperature. For windspeed, measurements 
were dependent upon different technicues of calibretions 
according to the windforce. All measurements used in this. 
popecrawere twenty-=minute means. The reader is referred to 


(2 for more cetailsweconcéernmings the instrumentation. 


Oo. Packground Theory. 

Extensive use will be made of Prandtl's turbulence 
model, because of its relative simplicity. This model is 
described in many stancard texts. According to this theory, 


tae following major results emerges 


' OU. 

ys = hb SO (1) 
¢ 4 4 % au \2 
Sees - pie <fiiies) (2) 


where W ,fr', are the turbulent velocity fluctuations in 
the direction of the mean wind and in the vertical, respece 
pave. Y ie is the eddy mixing lensth, while Ke is the roote 
mean-sovare mixing lensth. The eddy stress T' is then 
considered to be very nearly a constant C, within a layer 
ealled thewsurface layer. 

In the case of the fully-rowek, neutral surface ia yer, 
Rossby has shown thaths kz, so that (2) may te intergrated 
accom a lye 


= Uva. 
fl, = St fn S 2) 
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where Uy,2VTe/P is the friction velocity in a neutral 
surface layer, Z is the rovehness parameter and k= 0.38 
is the Von Karman constant. 

In the non-neutral surface layer, numerous theoretical 


profiles have been derived. However Deacon has proposed the 


semi-empirical relationship 


Ow _ 
Oz su (4 ) (4) 


Meere the friction velocityiUg.aNT,/P is constant in the 
non-néutral surface layer. The paraneter & is a decreasing 
munecion of Richardson umber. seictually Davidson and 

Tarad [4] have shown that 16 -1\ increases somewhat with 
heicht. However (4) still atfords a useful tool as an 
overall or Gulk=-relationsbip,. We useful resv1ts may be 
derived from (4). These are; 


(i) the intersrated wind profile: 


| i~f 
= gata [t&) - | 7 


(ii) the value My, by the limitin- process indicated 
in equation (7) below 


&u/5 (hz) 


My = k dim OS 
: EZ 7 Z, C2720)" 


so that 


KAS ee L Kim, So (7) 


Ze, 5 $a) 


Actually the lowest level at which wind data were available 
was Z= 25 cm and Z is the roughness ja rameter determined 
under the time-nearest neutral conditions. u, may therefore 
be approximated ty the finite-difference expression 
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joe = 0.38 Mas 


Im 25 - Wm Zc () 


Following a procedure similar to that employed by David- 


son end Farad [4] in neutral cases, it may be shown that 


Noowtt = aiee a(tnz) 


6 7 64u/F w; (9) 
where ee 
= Sco 4 
A(lnz) = In S92 + In 492 + Dy 209 (10) 
and 


ALL 


(ug-u,) #(us- Uy) 4 (u,-3) (ano 


for mean-wind measurements at 25, 50, 100, 200, 406, 800 cms. 


nee, GCeostromnic Drag Coemiacirentsr 

Lettau [7) “as defumegethe Meostrophic drar coef iicren. 
as a 
C = iye/Ny (a5) 
For neutral stability, the corresponding drag coefficient 
is defined as 

Coe Mo /Vo (13) 

Neutral wind, profibesmanme'Neill may ce obtained by findimc 
those cases for which the potential temperature profiles 
we characterized by 00e/0z=0 in the layer 25 cm - 2CO cm, 
Altogether 15 such profiles were found. Assuming a loearith- 
mic wind profile can be fitted to the wind data, U,. of 


equation (3) may be obtained as 


kaw — kak 
na SSO Ghee e ls 
o 





By analogy with flow in circular conduits, Lettau [+7 
has suggested a drar-coefficitent Cg for neutral csses of 


the form 


Ci gee Soe OHO 15 
Log (CaRe)= 2.24 wae) 


where no may be called the "surface Rossby number" defined 
as the non-dimensional expression 
Re =/fz , = Sh (16) 

The constants of equation (15) were obtained by Lettau by 
least-square methods applied to a wide variety of neutral- 
wld pYrotmiekhegeatlom diffemgent simttas, 

Lettau [7] goes on to give values of Cie for different 
values rtiyo9, the Richardson number at 100 cm, for the 195¢ 


O'Neill data. He obtains a curve as shown in fig. od. 
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Ratio of geostrophic drag coefficient for non-adiabatic and 
adiabatie conditions as a function of Richardson number at 
Z=100 em, at O'Neill, Nebraska 1953 (After Lettau [7] ). 

It must be pointed out that C, here represents the geostro- 
phic drag coefficient for any surface layer of drag coeffi- 


Crent C, Ll this layer were “suddenly tGovre convert d 2100 
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a neutral layer wi thout changing Ve/f Zo: 

A primary objective of this paper was to test eouation 
(15) as a means of determining Cg based on the 1956 O'Neill 
data. Also application of non-neutral vorofile data to eouae= 
tions (7) and (12) afforded values of C under identical 
newemmess eonditions. The comparison of Lettam"s* 1953 *g@rarh 
og C/Cy versus kKi1090 with that obtained by the writer at 
RLy49 (to be more exact, the Richardson number at 141.4cm) 
pawied on the 1956. 0'Neiiae data ismshewn in aig mss 

If one knows a value of Ro = Ve/fz,, it is possible to 
solve equation (15) for C, by numerical methods. It is 
shown in the Appendix that a numerical solution Cy, of eoua= 


tion (15) may be obtained iteratively from the eouation 


| y; 9.3026 t - + Ye + 6.5927 —- 2.3025 leq R ) 
Y,,. = ¥; a OE (eZ) 


where yyz i-th iterant to 0.2295/0,. The first estimate 


for Cy, was obtained from Lettau's curve [7 , Pigs 2. 
o.O30 


Ca. 


O.C4s , 


C.8024 


C. = ——_— 
flag (CR) - 2224 


0.04a 


e.C39 
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Geostronohic drag coefficient under neutral condition as qa 


function of surface Rossby number (After Lettau [7] ). 
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The chief inaccuracy in this approach appears to be uncer- 
taimby* in. theimnumeritcalaconstancomappecmumemain. coquation(15), 
especially at low geostrophic wind speeds, Another source 
Cf error iss themmecess liven Ssomcercms cc Ceuta mmo lating 
Vowpetween™ two maps 6 ormic hours Spartvin order to get the 
estimated value of Vg at the time of micrometeorological 
data-run. 

A second approach to determine the relationship between 
C and C, was used. Based on 15 neutral cases at O'Neill, 
Nebraska, values of Cy were computed using eouations (13) 
and (14). Values of Cy, were then plotted versus log Ro, 
end the resulting best-fitting: curve compared to Lettau's, 
ungich is replotted onmthe same eraph. Themconparisonmwill 
menatscussed in later wsection) 

The Richardson number has been used very extensively 
as a measure of stability in the micro-=meteorolorys its 


definition is 
g 20/02 


Re = @ (du /dz) 


This number is proportional to the ratio of work done arainst 


(18) 


static stability to the work done by eddy stress. 

For instruments arranged at @€cual lorarit*mic spa cima, 
as in the 1953 and 1956 Q'Neill data, a convenient finite- 
aitieremee ferm of Rl at 42 27) 5 has been proposed by 


Lettau [5] 
2 G - 6, 
RG,)2 Safe as 
2 ’ 


This value is understood to be applicable at the geometric 


mean of Z, and Zo, namely Z1.2=NZ12Z0. Qo and Lo in the 
8 ; 
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1956 data are the 20-minute mean ynotential temperature and 
windspeed at level ae an analogous meaning applies for 
symbols with subscript one, 
Finally, Monin and Obukhov [2\have introduced a use- 
tt mon-dimensionel paraneter Seder ime by 
G _ kz oi /2z_ (20) 
Uy, 
ieee tLnat the numée@ator of@s as the form 
kz a1/dz = Lp ai /dz (21) 
which, according to the Prandtl model, may be regarded as 
mepresentings Use. Une as here defined represents that 
part of Uy, due to turbulent fluctuations associated with 
tag mecrranical turbulent Wiximy Length Leake. The expression 


fore may then be written 


_ Map /Ve Ge 292 
>= Ux /VM ~ C me? 


ier cras coefiitient defined by Cr = Une/ Vg iswnoe ome 
that can be measured directly, tut is one whose value can be 
inferred ty empirical information re@arding S. This is 
dealt with in more detail in section 8, This interpretation 
Of Uxe as kzou/oz appears to have been riven first by 


Businger [3]. 


5, Comparison of C/Cqg for 1952 and 1956 data. 

In computine 0,, Lettau's empirical formula, ecuation 
(17) was first used with the 1955 data. It was necessary 
to compute Zo, using eouation (9) from the time-nearest 
neutral wind profiles and Vg from interpolated surface 
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cherts. In comruting Us,and C, eouation (8) and (12) were 
used. The value of Rij4o9 is also included using equation 
(19). The results for all data-days employed are given in 
me@ole 1, and the @foupings of Riqzo im classes are given in 
table 2. The graph of C/Ca versus Rij4o is given in fig. 3, 
with Lettau's curve \7) superimposed. The agmeement with 
Lettau's results is not as good as might have been expected. 
For instance the value C/C,g= 1.75 was found at Riz4o=0 
contrasted with Lettau's value C/C,=1.00 at Rizqo=0. More- 
over the maximum velve of C/Cg found in this paper exceeds 
Lettau's maximum by a comparable factor but with R1lj40 = 
-0.055, compared with Rly909= -0.007, according to Lettau. 

Possible sowrcest of error Ife in thie values of We, 
which affect the values of C and Cy, given by equations (12) 
and (15) respectively, as well as in the value of Z. The 
values of Z, determined here ranged from 0.25cm to 2.1 cm, 
as contrasted with values quoted by Blackadar et al [2] of 
Zo= 0.6 cm for profiles in approximately the same data- 


period. 
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PABLE I - Values of C/C, by equations (8), (12), (15), 
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Date CST Zo (ems) Uyem/sec) Ve (kts) C Gs C/Ca Rizao 
ema] I55 
LZOS eae 24,67 VESO° Omecs, Geese vos -0.46 
1405 vee 09.285 12.30 Gross Gree. Vie -0,17 
1505 1. Ee 01.24 12.00 Gl0ee Ones lee -0728 
LSS6O Lr, Le cOm9 ie soe 04040 Ga ne 120 @O,20 
1605 eke cOeld 11.60 0.056 O8024 1,54 -0O.20 
1705 fee cago 10.80 O7065 Giese Wa? -9.0 
1305 ie, Le 44,00 We52 0.084 Giese Zeal -C.C2 
els ee 2 lee 1z.CG Ose0E5 OF 064 1.02 0,06 
poroul 1956 
0950 o.14 49.29 7744 0.001  O, 087) ieee ea lo 
tio 3 oo14 42.835 1900 O1O8#5 Oye, Ieee! 
1205 Ores 46.69 Peo G2 OBZ On Gia? 1.939 -0,15 
1505 oe14 Denee e000 O.0S2 Giese lees -cye7 
1505 On ee 5 eit PeueoO Om os O-O65 aes =-QO,14 
1705 Oo. 14 eel OS) 16700 O8075 Geese Tey -iee- 
1805 e214 Gece 19.80 0.085 0.058 2azo =ORe2 
cléO oo 14 eee 9.52 0.059 0.040 1.47 0.04 
2c05 oe 14 Cleoe teoe O2085 0,042 “2Ael VCrGe 
Bo dul 1955 
0205 1.49 2215 £650 0,009 0,040 Oie22 Ome 
0405 eae Lens eeu OeOco —@.05/ O42 seme 
0805 1.49 90.49 24200 0,042 0.004 1teae —Ope 
O905 eo 52a Zo. OC GsO20 O20 ee eee -Q,01 
ie Onno 08 . 87 wa,00eC.028 0.029 Tay [ome 
1405 OO 69.24 “6e0C 7 02028 01.028 Tez Sere 
Leos Oe Te 6¢, 80 S400 502s GYO27 Ome 5 a), O | 
Boo O70 48.72 SH SOTO. ce 0.028 O70 saemer 
6 Aug 1956 
1805 0.55 Coes Ons Os ©9596 0,060 lea aeree 
2105 0.39 LAPS ROSS PNOROce O,OCO Gli wome 
p00 0.55 20400 PrevOm Oe055 O0.O0S0 TyoCsRCee” 
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6:7? Wevised formula for C, 

Hitherto Cg has been computed by using Lettau's empirical 
function for'C,, equation (15). Since Lettau himself reslized 
that the numerical constants in (15) were not very exact 7] 
an attempt is made here to obtain an improved relationship 
berween C, and Ro by using only 1996 O'Neill data. The first 
step again consists in computing Z, for the time-nearest 
reutral wind profile according to equation (9). As before, 
Use results from equation (14) with k=0.38, so that 

= Coot AL (25) 
Eouation (13) is then applied to determine Cg. The results 
Grete computations are shown in taple 36, 

An empirical function for C, was then obtained by the 
follrowin= procedure: 

Let log ro=zx and C,=y and suppose the desired Cuneo has 
cme f£ Orn 

y =be” (24) 
where y means the desired function, while b and € are cone 
stants to be determined. Takings the log of (24) one obtains 

Ye y = Qu b +e XK (25) 

Momswmolify the notations gen Yn y= Y and Lwhec , and there= 
fore the desired function becomes 

Y=ezar+cr (26) 
Using the least-squares method, the normal ecuations are 


se 


- Y-a-ck=0 
~ Cem 
D Yix; - a.) x; - cDx20 
b2. 24 tsi 
Substituting the vaiues of observations (table 3) in the 
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normal equations one obtains 
ao, Ole peor oe acl = © 

~58190-103..28a = 785.18¢c=0 
irom Which, it follows that 

@= -0,89 er Dae, 

Cee. 2570 
Substituting these constants in (24), and further by chanr- 
maeey CO C,, x to log Ho, the expression for Cy, men ve 
obtained 
Ca = @-427 Exp fe 0.369 heg Ro | (23) 

The graph of equation (28) is shown in fis. 4. compared with 
moubam S praph. From fig, 4 it Ws evident that the rewmisca 
formula for GC, tends to underestimate the drap at lary eq 
values of Rossby number, and overestimate it at smaller io 
(log to< 7.1). It should be noted, however, that there was 
considerable scatter of the observations relative to the line 
o. best fit. 

If one assumes that ecuation (28) is a universal formula 
form@mreostrophic drar coefficient itn widdle and hayher Watitude- 
under neutral conditions, one can then write 

Uae ait nz = SLE (29) 
with C, riven by (28). Hence, adiabatic wind speed at level 


Bis obtainable as a funcUion Of Vel Z4, 2 anda lace ge ee 
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Geostrophic drag coefficient under neutral conditions as 
a function of surface Rossby number. 
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7. A revised formula for C/C, 

Themrevised foemmula forse, wil) besusedm Geis section 
to compute C/Cyg, still vsinse the previously tatulated values 
of seostrophic winds and C (table 3). The results are shown 
in table 4 snd tre resultins values cf C/C, then plotted 
wersus Rijyay as shown in fags o. The procedure of fitting 
inpe Curve is shown below. 

After plotting the values of C/C,, a furtrer restriction 
is placed on the empirical function C/Cg: namely that C/C,=1 
at RicO. horeover, since the shapes of the curves for C/C, 
in unstable and stable conditions are different (according 
to the scatter cf points), it is assumed that in stable condi- 
tions the desired function has the form y= e?* where y=C/C, 
and x= itljy4o- In unstable cordition, y is assumed to te of 
the form y= \+ k we” where y=0/C, ane X= =i AG, sae 
ko, are constants to te deter:ined. 

(1) For stable conditions (table 5), with Riyao= x and 
O7ea= y, if the assumed funcivoneaaeectoeicrm 

v= e** (20) 


then 


end if Inn y = Y 


By the method of least-scuares, th§ normal equation will te 


Vv Mm . 
dD. Nix, -a2 et (2) 


- 6 
iz) t=3 


Solving for a, 





Sukbsbituting this constant in (30) and changing y to C/C, 


BLO o RL 405 the sr oeesslon Lop C/Ce ls obtained as 


C/Czme Cn . 4.159 Rissel (z 





TABLE 4. - Computations of C/Ca using equation (28) for 


Date CST 


10 July 1956 


co 


25 


1305 
1,05 
L505 
1530 
1605 
1705 
1805 
2105 


July 
0930 
1105 
1205 
1305 
1605 
1705 
1805 
2130 
2305 


July 
0205 
0,05 
0805 
0905 
1230 
105 
1605 
2330 


1956 


1956 


6 Aug 1956 


1805 
2105 
2205 


7 Aug 1956 


O05 
0805 
0905 
1105 
1330 
1605 
1805 
2005 
2205 


log Ro - 0.369 fog Ry e 


6.799 
6.773 
6.7h5 
6.&60 
6.730 
6.699 
bon / 
6.75 


6.500 
Cee 
6 .hi6h 
6.139 
Geoc 
6.1113 
6.107 
Bo 30 
6.084 


6.215 
6.516 
ep fel 
6.948 
7.366 
7.523 
7.593 
e512 


7.18E 
71.96 
7.20h 


6.295 
6.330 
CaS | 
6.350 
Peat: 
TeOe? 
Taloc 
fae ea 
7.000 


-2.51 
-2.50 
-2 9 
~2.53 
=2eine 
~2,h7 
-2 16 
~2.)19 


-2.110 
-2.39 


eo 638 


=e 
-2.36 
2.36 
-2, 36 
— 20 
eZece 


~2 29 
=2 110 
2,17 
-2.56 
-2./2 
-2.77 
-2.€0 
-2.77 


~2.65 
=OO0 
-2 66 


~2.32 
“2.33 
-2.34 
-2.3 
-2.62 
=D Ol. 
-2.65 
-2.60 
-2 56 


p 


0.081 
0.082 
0.083 
0.080 
0.08) 
0.085 
02065 
0.083 


0.091 
0.092 
0.093 
0.093 
0.09) 
0.09k 
0.09k 
0.100 
0.109 


QcdlOI 
0.091 
0.0€&5 
0.077 
0.066 
0.063 
0.061 
0.063 


0.071 
0.070 
0.070 


0.098 
0.097 
0.096 
0.096 
O<One 
0.07) 
CO. 
0.07) 
Oneal 


Ca 


0.0348 
0.0350 
0.0360 
0.03h2 
0.0359 
0.0363 
0.0363 
0.0360 


O2@305 
0.0393 
0.0397 
O20377 
0.001 
0.0Qh01 
0.0401 
0.027 
0.065 


0.031 
0.0388 
0.0363 
0.0329 
0.0282 
0.0269 
0.0260 
0.0269 


0.0303 
030277 
0.0299 


ole 
-O41S 
2010 
FOO 
soe 
00316 
.0306 
20316 
O.0Se27 


eel aL eC) (57) © 


ao 


C 


0.0363 
0.0560 
0.0521 
0.0396 
0.0530 
0.0686 
0.0838 
0.0350 


0.0507 
0.051 
0.0515 
0.0623 
0.069 
0.076 
0.0853 
0.0591 
0.0&50 


0.0090 
0.026) 
OcOu2d 
0.002 
0.0368 
0.0284 
7.0233 
ences 


0.0557 
0.0216 
0.05) 


0.0970 
0.0979 
0.1666 
Oeiesa2 
0.0755 
0.06 64 
0.0437 
0.031 
0.0153 


C/Ca 


1.05 
1.600 
Lelbs 
1 rGo 
1Lel75 
1.785 
2.320 
0.975 


1.305 
1.145 
i OS 
1570 
1.615 
1.860 
2.150 
1. 3R5 
1.830 


0.209 
0.680 
1.60 
1.220 
1.305 
1.055 
0.895 
0.810 


1.835 
0.722 
1.515 


Cane 
223600 
1.060 
2.510 
2 li20 
2 co 
1.540 
1.069 
0.65 


Ca 


Rij Lo 


~0 6 
~O.17 
-0.28 
-0.25 
-0.20 
2s 
=0.02 
0.06 


=“0015 
-0.21 
=0.16 
-0.09 
-O.1h 
-0.06 
-0.02 

0.0h 

0.035 


Gal 

0.05 
=0.10 
-0.01 
=0.02 
~0.0? 
=(@\, 61. 
0.00 


eO00 
07 
Ol 


Ore oS 


0.01 





TABLE lh. (Cont'd) 


P= 
Date CST log Ro ~ 0. 36ghesR e? Cc C C/Ca Rin). 
27 Aug 1956 
1230 70390 <=2.73 (02055 “G.0275 Og@e25. 15976 5-0.07 
1505 7.361 €©2:72°°0.066 0.02829 070307 M070 = —0.02 
1705 7.308  =-2.72 .. 0.066 0.0282 0,0316 17975. ~-0.02 
1930 7.168 =<2.64 0.071 6.0308 088023 07676 0.15 
2205 7.392 <-2,12 "@.066 “0.0282 O50022 SO -ene 0.0) 
29 Aug 1956 
2135 6.377 -2.35 0.095 0.006 0.0669 1.60 0.03 


TABLE 5. - C/Ca in stable conditions 


L Rage k: = CKCazy = Imy, 2 Ys YE x; th 
ll 0.057 0.995 -0.925 -0.001 0.0C3 
2 0.036 1. 386 0.326 Geol 0.001 
3. 0.035 1.830 0.60) OyOen 0.001 
he = 0069 0.209 -1.565 =0.155 0.010 
oO rO5 0.680 -0. 386 =0.020 0.003 
6. 0.067 Oia -0. 326 -0.022 0.00) 
foe, OLOl! ioe 0.415 0.006 0.0CO 
cy 07000 Pes els 0.2 G.007 0.000 
GS. 0.05) 1.069 0.068 0.00) 0.003 
HG. O-127 0.165 -0. 766 -0.010 0.016 
ie. yO. 163 0.076 -0.27h -0 .0))2 0.023 
eye 0ieOh0 0.078 -0.2h8 -0.010 0.002 
ie. 0.031 1.60 0.495 0.015 0.001 
Suns -0.283 0.068 


ot 





= ee se 
aii teen. 


> 
- ete t ‘ 


(44) For unstable conditions (table 6), the scatter 
of points of C/Cg versus Rij49 suggests a curve of the 
form 


yei+ kx e%* (22) 
where - Rijzc=x and C/Ca= Yo 


Again C/C, has been forced 
to approach unity as RL approaches zero, 


If we let y-l=¥ 
and then take natural logs 


Bmp = Unk, + Inn + hye (34) 
Jn, © = Mek, + kv 


Let Ind = Y and Ww = a. , and substitute in (34) 
ve = (oe ke 
The normal eouations are 


Y 2.0 BR xc 


NW iM A , 
DY-x: = a.) Xi - ky 24 =o 
ts4 isa : 


(35) 
Solving for kj and Ko 
Kats oO, 564 9 Ko = e135. 4204 
Substituting these constants in (23) and replacing x ty 
mriac, yy by C/Cg, gives 
GiC,. et. — 30.364). exp [13.324 Rigo] (36) 


Note that C/C, of equation (36) may be maximized for Riyz,4¢ 
= -0,075, for Which (C/Cz)),. is 1.86@. Mats Mie 
close te the nfraximum of the fitted curve of fig, 5. Letvau's 
graph of C7C 3 has been included in fig. 5 and it may be noted 


that the computations of t'ils teper are in general agreement 


fe 





with those of Lettau, although his data did not span as large 
a rance of Ri values. The results indicate that C/C,)1 in une 


stavle conditions, whereas C/Cg <1 in stable conditions. 
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TABLE 6 - C/C, in Unstable Conditions 


Y=Clen Pez yor 


1.05 
1.600 
1.5 
1.160 
1.475 
eae 
aaee 
1.305 
1.1h5 
1.305 
1.570 
1.615 
1.€60 
2.150 
AGO 
1.220 
I sles 
1.055 
2.360 
2.510 
2.1120 
PO 
1.5h0 
70 
1.090 
aes 


0.0h5 
0.600 
OuhS 
0.160 
0.4.75 
0.785 
1.320 
6.305 
0.1)5 
0.305 
0.570 
0.615 
0.€60 
1.150 
0.160 
0.220 
0.305 
0.055 
1.360 
ASA 
ieee 
1730 
0.50 
O.1 76 
0.090 
0.125 


¥./x; 


0.097 
3.630 
1.590 
0.638 
24570 
14.700 
77.600 
2.080 
0.690 
1.90 
62.600 
l. 360 
14.800 
52.30 
1.550 
20.000 
17.950 
2.890 
32.400 
19.600 
33.000 
157. 300 
514.000 
22% 100 
3.750 
5.200 


wh .¥ 


=i9 253 


1.289 
0.6h 


-0..Nh9 


Ovo 
2.688 
4.352 
0.737 


= ee 


0.663 
137 
Le72 
2,695 
BR OST 
0.38 
2.996 
2. GCG 
Ie eew 
3.478 
2.976 
3.500 
5.056 
3.989 
3.122 
le 322 
1.69 


52.6h9 


x" 
Cae 
0.027 
0.07& 
0.063 
0.00 
0.003 


SS 


0.022 
0.0hh 
0.025 
0.008 
0.020 
0.003 


0.011 


Vix: 


=] © O80 


Cont? 
0.130 


“0.113 


Ov 173 
0.12 
0.07) 
0.108 


-0.078 


0.10); 
0.376 
0.208 
0.156 
OS 007 
0.045 
0.033 
0.09 
0.020 
0.16 
0.229 
opisal 
0.056 
0.00 
0.23h 
0.032 
0.00 
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Ratio C/C, as a function of Richardson number at Z=140 cm, 
at O'Neill, Nebraska 1956 with Cy, given by equation (28). 
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8. The normalized—lorarithmic wind sheaz S 


It was explained in section 4 that if one knows 


L- hz du/dz 
—— 


Ellison [2] has given a prediction equation for S in unstable 


one may obtain directly Cr/C,or Uxe/Wee- 


cases 


' a , 
s+ T25 = | (27) 


Wimewen , 1s helLeht and L is the @radient length defined as 


pia _ Une OU Soe (36) 
kg 08 /oz 
with YY cGefined as 
eauela | 
= _—_ 349 
x K.. y (39) 


here © is a constant, and K, and K,, are eddy diffusivity 
coefficients for heat and momentum respectively. Dividing 


toth sides of (37) ty 34 


-& 
Substituting (20) and (18) in the expression for a and 
Ls 


Solving for 8 


ee 
S = (1-78) _— 
Blackadar et al [2h suggest the value ¥'=18 for urstelle 
Corncitions. However in stablewcendiiiens, they fimausea 
equation (40) gives no tetter approximation: than the hMonin- 


Obukhov "log + linear" wind profile defined ty 


S — | + ae Zz ; X =—CoOnstant 


which is ordinarily valid for Richardson numbers near one, 


(41) 


Apain, it may be verified that (41) is eouivelent to an 


ao 





equation 

Se itn (42) 
Blackedar et al [2] have tabuleted values of &’ at C'leill 
fOr Varlouseiinvers. Por the laveraiecureseer stable con- 
Gitions, the best=fittinge value om the Cemaweant was X's + 
which widll=be used for the stam@le cascs Sings) paper. 

The computations for 5 are Saown in tatile 7, The values 
of S are presumed applicable st the heicht assigned to Ri, 
peat 15 ateZec V2 m or Z= 41.4 em, The curve ome 
shown in fig. 6. The result indicates that S¢1l in unstable 
Gases, and S>1 in stakle casese Whether S&S in stable condi- 
tions levels off with a further increase of Ri, is not cer] 
tain yet, since the value chosen for &' may pot e valid for 


are DOSitive values of Ri. 


a 





TABLiG./. = Results, of conputat ons formes. 


Date CST RL 140 Ve (kts)  $(y40) 
1Oeasly 1956 
1305  =0.46 z.60 O'eevie 
1405 -0.17 12.80 0.708 
1505 -C,28 12.00 0.628 
150 -0,25 15.64 0.652 
1605 -0.20 11.60 0,683 
1705 -0.05 10.80 0.850 
1805 -0,02 1) 32 OnGe5 
2105 0.06 12.00 Vy SLS 
23 July 1956 
O9Z0 -0.15 19.44 0.724 
1105 <Oret 19.00 0.676 
1205 -0.16 18.00 Ces 
1305 -0.09 Lee 0 784 
1605 -0.14 16%20 0.729 
1705 -0.06 16.00 0.826 
1805 -0,02 15.80 0.916 
2150 0.04 9.52 lee 
2305 One Ta 52 Teele 
25 July 1956 
0205 0.10 4.80 1.888 
0405 0.05 9.60 ay 
0805 -0.10 24.00 0.769 
0905 -0.01 26.00 0.958 
1230 -0.02 32,00 0.936 
1405 -0.02 46.00 0,929 
1605 -0.01 54.00 0,940 
= 0.00 44.680 AOSS 
emnue, 1SOG 
1805 0.00 OED 0,994 
2105 0.07 10.80 1.805 
2205 0.01 11.00 ieee 
Veanac 1956 
0405 0.01 11.88 1.076 
0805 -0,04 12.88 0.868 
0905 -0.02 Tess 0,929 
ie. -0.08 15,48 0.804 
ie a@ -0.04 15Ke4 0.867 
1605 -0.01 1x00 0,958 
1605 =Ongil eo 0.940 
2005 0.05 14.00 1.420 
2205 Ole 12.40 2.147 


ee 
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TABLE 7. 


R140 


-0.07 
-0.02 
-0.02 
O.19 
0,04 


(Cont "ap 


ao 


Ve (k t 8) 


26.48 
24.80 
20, 20 
ec 
26, 60 


5(140) 


QO, 808 
Cg 
0.914 
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big. sof 
The normalized logarithmic Wind shear S as a function of 


Richardson number at Z= 140 cm based on 1956 O'Neill, 
Nebraska data. 
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9, An alternative method of cetermination of Uw 

in Securtn 7; c/C., was computed by eovations (8), (12), 
and (28). While there may te considerable error in assuming 
a funetreonal Tom for@™¢, ofthe type dermucomea) Securen >, 
Enere is alse a possibility of error in tme=ecenouration of 
Ux by equation (8). Hence an alternative computation of Uy 
wes employed, one which employs the expressions ior 3 of the 
section 8. This approach essentially leads to a value of ‘A, 
in terms of the wind-profile UC’) and the intesral-mean value 
of S in a suk-layer of the surface layer. Thus for example 


S a hz OUL k Ou 


“Uy O2 Uy Sinz 
and integrating from Z=Zo to Z=Z, assuming Uy, is constant 


in tke surface laver leads to 


_ Us | Zz 
a 5 Yn. Z (42) 
© 
where S$ is the interral-mean value of S up to level Z. 
As in the previous section S was taken in accordance 
with 
th 
S=(i- 18 Ret2)) unstable cases 
(44) 
S= '+9 Rice) stable cases 
elenowah Priestly (8] has raised a question regarding the 
the validity of a constant Uy, with heirht in stable condi- 
tions, it was assumed that a surface layer at the depth of 
at least 2m exists for toth stable and unstatle cases, 
Values of S according to ecuation (44) were plotted 
using iti (Z) for velues of Z= 5, €5, 7O, 140, 230c. aac 


lowest value of Freight corresponds to a level at the ®eorietric 


Cn 


1 





mesn of the roughness parameter Zo and Z=c5s cm, weime the 
cnown values of OU/d2 and 36/dz for this layer. The 
mraonis of 5 versus kn z was plotted and the intestral-mean 
value S extracted. Fip. 7 shows an example of an unstable 
end stable case which occurred at 1905 and 2lOs CsT, respec- 
Gively on 10 July 1956. This permits one to solve for Uy, by 


means of 


k U(200) 
= Yn (200 S20) 


The values of Uy computed by equation (45) are displayed in 


—— (45) 
table 8 alonp with the values of Wg, listed earlier in 
table 1. The absolute averare error is approximately 10%, 
althouch the algebraic error is 3%. 


Note that the ratio Cr/C,g may te obtained as the product 


cS. = 
Ca. 


The value of S as a function of Ri(Z) may be ottained from 
equation (44). Values of Uy may be obtained from equation 
(45); and C, may te determined from eouaticn (28). 

The main purpose of this last discussion is to empha- 
size that Cr is not equal to Cg. Listed below in takle 8 


are some computed values of Cree 
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The normalized logarithmic wind shear S as a function of 
Lenz for unsteble and stable conditions for 10 July 1956. 
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TABLE 8, Cr/C, in unstable and stable cases 

Date CST Ri14o s Ue: /3ec) (Ce/Cg), Ugnsec) from (3) 
10 July 1996 

S08 -O0.46 Ceo Cone OCs 24.7 

L405 Sere ee 0. 9CO 59.0 eS WO o oI 

005 -0.28 O.E85 Ole S 0.925 oie 

605 -0,20 O.9CC lee 1,028 ono 

Lies -O.05 Opeoo 0458 1; DS cee OC 

La05 -0.C02 Oe 42,0 ea OU 44,0 

ZrOD O06 ee SO ZU Loaes eo 
23 July 1956 

OS -Q.21 Dee oO ote mi 0,648 42.9 

205 -O.16 0.880 Zoe! Ogres 46,7 

ee -0,09 OFES5 Cree Povo 52 eee 

EOS -O0.24 Omseds 4205 OO. 956 5 eS 

17S -0.06 05750 45.7 belo oes 

Be09 0.04 Poo Fo Cae. cleo 
25 July 1955 

0405 Oreos ie POO ion ee Wah 

Aus 1956 

2105 O07 Nines ZCO Ley) Lee To ey 

Aug 1956 

ZOOS 0.05 POO eto Poot ere 
27 Aur 1955 

Dae 0.04 1.00 FOS OncuG Cae 





10., Summary and conclusions. 

A functional relationship between C/oe and Ri was 
edotained by SWsime@ean expression, for™eg™surcccoveusoy Lettan, 
Better agreement was obtained after making some revisions 
in Lettau's formula for Cg. An expression for C, has the 
oractical significance that one™@¢an.Bet directly the ratio 
of a windspeed at some anemometer height to a surface feo- 
strophic wind by knowing the roushness parameter. 

By using the revised expression for Ca"°4 new iumeerema. 
relationship between C/C, and ri was obtained. This function 
Ibes different analytical forms in unstable and stable condi 
TOMS « 

Another way of obtainines C/Cyg was based upon computa- 
tion of Ux using Blackadar's results. By this means a 
wind-profile expression was obtained which varies with 
Stability. From this new precedure Gene can met the valve 


Bimet/G, “itnout using C/On: 


OO 





APPENDIX 


Computation for Cg using Lettau's formule 


The original ecuation is 
0.102 
Ci Se ee 
Los (Ca Roy = 2.24 
Ey eeross multiplying 


Log (CaRe) rma 2.24 + O-104 TCE 
or | (2.24 +0:104/Ca) 
Ca Rese 
0.104/C,, 
= 173.8 x 10 
then Aes 
O.I1C 4 rw 
Cake — ci hl 
173.8 Ca 


2 Kjo ( 0.23g5/Ca)} 
and substituting in the above eauvation 


¥ a y 
Fd) = Re ( aa =398 | oo, =o 


(3 


Met Ose 595/Ca= y 





173,8 ¥y 
- 0,2395/173.2 Ro — ye” - oO 
7 y 


= Tace VV #0 


y 
. oo *, = 9 pe C- 
Fry) = 0.2395/173,8Ro - ye =O as 
Taking the lors to the base 10 





0.24399 ; re 
Ket tor (SSF R= kK —- constant, and replace log y by 


0.4242 In y and substitute in (48) 


Fry) = 04343 hy + 0.4343y-Kaog (49) 


To get the root of ecuation (49), we will use the 


On 


nm 
SD 





Newton- Raphson method 
Cy, 
yr = a s Bid (50) 
; , 
Since f (y= 0.4243/y+ 0.4242, therefore by substituting this 
in (50) we obtain 


Veep = Y, — 014343 bays + 0.4343 ys — IK 
O-4343/Y; 4+ ©-44343 
Substituting log (2772? Ro ) for K and divide the fraction 
by O.4c4e , 
0.42436 
ie si V,: = Yn q es as \/o.4 343 bea 
Jy + 





Changing Imy, to 2.3086 log yy 


mee 2, Cie bog a + 6.5927 — 2. 3015 Leg Re (5) 
ao aS ae a ene 1) 


Ol 
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